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Relativistic Linear Restoring Force 

D. Clark, J. Franklin, and N. Manif] 

Department of Physics, Reed College, Portland, Oregon 97202, USA 

Abstract 

We consider two different forms for a relativistic version of a linear restoring force. The pair 
comes from taking Hooke's law to be the force appearing on the right of the relativistic expressions: 
^ or ^ . Either formulation recovers Hooke's law in the non-relativistic limit. In addition to 
these two forces, we introduce a form of retardation appropriate for the description of a linear (in 
displacement) force arising from the interaction of a pair of particles with a relativistic field. The 
procedure is akin to replacing Coulomb's law in E&M with a retarded form (the first correction in 



^ the full relativistic case). This retardation leads to the expected oscillation, but with amplitude 

Ph 

growth in both its relativistic and non-relativistic incarnations. 
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I. INTRODUCTION 



By the end of a student's undergraduate career, the classic harmonic oscillator has ap- 
peared in almost every imaginable context. We see it in introductory physics, as our first 
example of a spatially varying force, again in a sophomore-level physics class, as a vehicle 
for studying the Fourier Transform, and maybe again in a computational physics course, as 
a model for bonds connecting atoms. Aside from its direct utility, the equation of motion 
associated with the harmonic oscillator occurs in a variety of settings, and most instructors 
refer us back to our earlier experience with masses connected by springs as a way to visualize 
widely disparate physical systems that have, at their base, the same fundamental equation 
and solution. One reason for the central role of the harmonic force is its near-universal 
applicability - for almost any potential that has a minimum, motion in the vicinity of that 
minimum is governed by a linear restoring force with a "spring constant" related to the 
curvature of the potential at the minimum. 

When we move to a relativistic description of simple harmonic motion, that universality is 
lost. There is an a priori ambiguity in the interpretation of a classical force F in relativistic 
dynamics. It could either appear as a component of a three-force, so that ^ = F, or it could 
be a component of a four-force, in which case ^ = F. In E&M, for example, the Lorentz 
force is an "ordinary" three-force (see p. 519 |IJ). For our non-relativistic F = —kx, we must 
first determine what type of formulation (three or four-force) to use. However, both of these 
formulations share the same behaviors in the non-relativistic and ultra-relativistic limits. 
Indeed, we could argue that correct behavior in these limits is the minimal requirement for 
a system to be considered a "relativistic spring." In fact, we can define an infinite family 
of equations of motion with the right limiting behaviors to be associated with a relativistic 
form of Hooke's law. Following [2J, we will focus on the three/four-force pair associated with 
F = —kx. We'll start by considering the relativistic dynamics of each of these. 

If we want either form of the force to correspond, as in E&M, to an interaction of particles 
with fields, we must introduce a retarded time evaluation. We proceed to do this for a pair 
of particles, assuming that the governing field moves at speed c in vacuum. We assume the 
force depends linearly on the distance between the particles, taking into account the time of 
flight for the field linking the particles. This process represents a minimal substitution, in 
which we effectively replace the instantaneous evaluation of displacement with a retarded- 
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time evaluation. For this modified form we find a physically implausible growth in oscillatory 
amplitude. The pair of particles by themselves no longer constitute an isolated system, and 
unless we include contributions from the field, we appear to have lost conservation of energy. 

This paper is meant to demonstrate the logic and physical checks that are common in 
high-energy and gravitational model building as applied to a toy system that is familiar 
from the undergraduate curriculum. In the spirit of demonstration, we take a pair of simple 
choices (letting —kx play the role of a component of a 3 or 4- force) and explore the physical 
implications of these two systems. 



II. NON-RELATIVISTIC SPRINGS 

Let's start by reviewing the standard problem and its solution. For a mass m attached 
to a spring with spring constant k and equilibrium at zero, Hooke's law combined with 
Newton's second law give the equation of motion: 

— [mx] = —kx. (1) 
at 

If we start with a;(0) = a and x{0) — 0, then the solution is 

<t) = « cos ?j . (2) 

One immediate (special relativistic) problem with this solution is that the speed of the mass, 

x(t) — a \ — sin ( \ — t 1 , (3) 
V m V V m 



has a maximum value of u^ax = \l — . This value can be made greater than c for some 
choice of a. 



III. RELATIVISTIC LINEAR DISPLACEMENT FORCE 

Suppose we start by considering the simplest possible equation of motion, given by ^ = 
—kx. We'll take p to be the relativistic momentum, giving 

d 
di 



mx 



1 - K 



—kx. 



(4) 
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From this form we can immediately make a qualitative observation about the extreme rel- 
ativistic oscillatory behavior. Assume the same initial conditions as above. For a large, 
the force is large, and the mass undergoes rapid acceleration from rest. But the maximum 
achievable speed is c (this is enforced by the form of the relativistic momentum appearing 
in Q), so the acceleration must die off as the mass approaches speed c - the larger the 
initial extension, the faster the value of c is attained. Once the mass has travelled to —a, it 
is at rest, and the process begins again in the opposite direction. Taken to the extreme, the 
resulting motion looks like a sawtooth, with the mass moving at a constant speed c from a 
to —a, instantaneously reversing direction, and again moving at speed c from —a to a. A 
physical manifestation of this motion would be light bouncing back and forth between two 
mirrors. 

In fact, any sensible "relativistic version" of a spring must have this behavior in the 
extreme relativistic limit, in addition to behaving like a traditional non- relativistic spring 
in the opposite limit. Aside from details, we have a clear picture of what must happen: 
a smooth transition from the non-relativistic sinusoidal motion to the extreme sawtooth 
motion. This is already an interesting prediction, since we know that in the non-relativistic 
limit the period of the motion is given by T = 2 vr ^/^■, a quantity that is independent of the 
mass's initial conditions. In the extreme relativistic limit (where a is large), the period of 
the motion is T = 4a/c, and so is governed entirely by the initial extension, with no spring 
information involved. One interesting consideration is the manner (as a function of a) in 
which the system "loses" its knowledge of the spring. We will find that the two forms for 
our equation of motion, ^ = —kx and 4^ = —kx, give different results. 



A. Hooke's Law Forces 

Consider a diametric hole drilled through a uniformly charged sphere. A charged particle 
moving in the hole experiences a linear restoring force. The equation of motion is: 



dp d 



mx 



—k X, 



(5) 



dt dt 
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where k is fixed by E&M constants. We can isolate acceleration on the left, giving us an 
equation that can be compared with its Newtonian counterpart: 

3/2 



mx = —kx I 1 — 



X 



(6) 



Starting from rest with extension a it is possible to solve for x{t) in terms of incomplete 
elliptic functions [3]. 

In we took —kx to be a component of the three-force, so that ^ = —kx. This 
formulation makes sense for a particle moving in an external harmonic potential, one in 
which the constant k is fixed in the "lab" frame. But, we could also let k remain fixed in 
the rest frame of the moving particle - in this case, ^ = —kx is the relevant equation of 
motion, with —kx appearing as a component of the four-force. We could engineer this case 
by changing the density p{t) in our uniformly charged sphere so as to achieve constant k in 
the instantaneous rest frame of the particle [? ]. Physical configuration aside, if we move the 
non-relativistic potential | kx"^ to the Lorentz scalar potential | kx'^ x^ (for Xfj_={ct , x)^), 
then the (spatial) equation of motion would read: 



dp 
dr 



dp 
dt 



d 
di 



The form equivalent to ^ is 



mx 



-kx ( 1 



X 



mx 



-kx. 



2\ 2 



(7) 



(8) 



To express both of these options in coordinate time, where ^ is the relevant quantity. 



we will write these two equations as ^ = Ft and * = F^, with Ft = —kx the three-force 



dt 



formulation and Fr 



-k X 



the four-force formulation. 



IV. QUANTITATIVE COMPARISON 

We now want to explore the details of our two dynamical systems. This is made clearer by 
rendering all of our equations dimensionless. To do so we use the length scale a set by initial 
conditions and the natural time scale ^/mjk (from Section |ll|. Let x = aq and t = ^Jmjks 
so that a;(0) = a becomes g(s = 0) = 1. We also define the ratio of initial potential energy 
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— k a 

to rest energy, a = ^^^2 ; a parameter that governs the size of the relativistic correction. 
Then our equations of motion become, in order (i.e. ([T|, (|6]), ([S])): 

?7v(5) = -(1n{s) 

q:is) = -qt{s) {l-2aq[{srf' (9) 
€{s) = -qAs) (l-2ag;(.)2)^ 
We can solve all three equations in ^ for various a - the result of that numerical 
calculation [? ] is shown in Figure [l} Referring to the figure we see that for small cr, the 




a = 5.0 a = 10.0 



FIG. 1. The non-relativistic motion is shown in black {q]\f{s)), and relativistic motion under the 
influence of Ft and Fj- are shown in gray {qt{s)) and dashed gray (g,-(s)) respectively. The motion 
is displayed for increasing energy (ratio) a. 

three separate equations of motion have very similar solutions. This confirms that both 
relativistic forms revert to the non-relativistic solution for small a, which is also clear from 
the expressions in (g. As a becomes large, both relativistic forms move towards a sawtooth, 
as they must. 

As expected, the low energy limits of the two trajectories coincide, and the high energy 
limits, even at a = 10, are beginning to converge. However, Figure [l] does not present a 
good way to explore the intermediate regime. The solutions to (|9| are periodic for all a, 
and we know the period of motion in both limits: the Newtonian limit 2 vr is immediately 
obvious, and the sawtooth period is 4 a. We can then compare Tt{a) and Tr{cr) (the 
periods of the motion under the influence of Ft and Fj- respectively) by direct integration. 
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The equations in (|9j) are all of the form 

q"{s) = -q {l-2aq'Y (10) 

for p = {0, 3/2, 2}. This general form can be integrated once to give 

^^(1-2^ 1 

4a(p-l) 2^ ^ ^ 

for constant A. In the case of p = 0, the constant A is the non-relativistic energy (with 



a constant offset), and for p = 3/2, it is the relativistic energy. For p = 2, (11) implies a 
conserved Hamiltonian, however A cannot be interpreted as the energy [? ]. If we use the 
initial conditions (^'(0) = 1, ^'(0) = 0, appropriate for a particle starting from rest) to fix 
the constant, then we can solve for q'{s) in terms of q{s), yielding 

1/2 

(1- \l + 2a(v-l) ll-a^]\'-^ 1 

^2 a 



g'(.) = ±-L (l-[l + 2a(p-l) (1-g^)]^)'^'. (12) 



Note that one can also consider (10)-(12) for other values of p as they define an infinite 
family of relevant effective forces. Any suitable relativistic model must have the same non- 
relativistic and ultra-relativistic limits (as these do), but beyond that, there is freedom. In 
addition to the two specific cases of interest to us, a model resulting in p = 1 was developed 
from material considerations in |'5|. 



Using the general expression ( 12 ), we can integrate over a quarter of a period with definite 
sign. For qt {jp = ?>/2), this gives: 



TM) =A^—[{2 + a) E{a/{2 + a)) - K{a/{2 + a))] (13) 

and for qr (with p = 2) we have: 

T,((t) = 4 Vl + 2aE(2(T/(l + 2(t)). (14) 

In both expressions, K and E are the complete elliptic integrals of the first and second 
kinds [? ]. These two periods, together with the limiting cases (the sinusoidal Newtonian 
one, and relativistic sawtooth), are plotted in Figure|2| We see that the period T^-iya) diverges 
from the non-relativistic period more quickly, and converges to the relativistic limit more 
slowly, than Tt(cr). Hence the transition between the two limits is more gradual for Tr{(y) 
than for Tt{cr). The period for any value of p should look similar to these, with a different 
transition from one limit to the other. 



7 



log 



2.5 




log(cr) 



FIG. 2. The period of motion as a function of a for the solutions to Q - here, we use \og{a) to 
capture a wide range of behavior. The non-relativistic and relativistic limits are shown in black, 
the periods under the influence of forces Ft and F-j- are shown in gray and dashed gray. 

V. INFORMATION SPEED 

There is another element to relativistic systems: the question of causality. Our model 
for exploring this will be a linear (in displacement) restoring force with a finite speed of 
propagation for the interaction. For concreteness, we will take the form F^-, although the 
qualitative features of the dynamics produced by Ft and Fj. are similar [2j [? ]• We will 
suppose that underlying the production of F^- is a field that mediates the interaction between 
two particles as in Figure [3j Given material considerations, it would be natural to imagine 
that field moving between the particles with some fundamental speed f , determined by the 
properties of the supporting medium. 




m m 



FIG. 3. A field mediates the interaction between two particles of equal mass, producing a force Ft-. 



If we have two equal masses that interact via F^-, and we do not consider a time-of-flight 
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correction, then the equations for the left and right masses should be: 



Fe{t) = k [xr{t) - xe{t)] Vl-^ 



Fr{t) = -k [Xrit) - X,it)] 



(15) 



where each mass experiences a force that depends on the instantaneous location of the other. 

In order to introduce dependence on an information propagation speed v, we will evaluate 
the force on the right hand particle using its location at time t, but use an earlier position 
for the left-hand mass: xe{tr), where tr is the retarded time defined by the condition: 

V {t - tr) = \Xrit) - Xe{tr)\. (16) 

Assuming the masses are given symmetric initial conditions (so that the left mass starts at 
Xi{0) = —a with xe{0) = 0, and the right mass has Xr{0) = a with 0:^(0) = 0), we can write 



Feit) = k [xriU) - xeit)] ^/l-^ 



Fr{t) = -k [Xr{t) - X,{tr)] \/ 1 - ^ 



(17) 



This move is analogous to taking the Coulomb force between two particles in one dimension 
(labelled left and right as above): 



Aneo \xr{t) - xe{t)\^' 



and replacing the difference Xr{t) — xe{t) with Xr{t) — xiitr) (with f = c in (16)). That 



replacement would yield the correct force up to corrections of order Xr/c that come from the 



actual field of a moving point particle. Motivated by simplicity, we will start with (17) and 
see what solutions emerge. 

If we begin with the stated (symmetric) initial conditions and define x{t) = x,.(t) = 
—Xe{t), then the single equation of motion together with the retarded time definition is 

mx{t) 



Pit) 



1 _ ^W! 

-L n 



dp(t) , , , , , f ±{{f (19) 

-^ = -k [X{t)+X{tr)] yl--^ 
V {t — tr) = \x(t) + x(tr)\- 
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In order to get the right-hand side of the equation for in terms of p, we must invert 



the first equation in ( 19 ) to find x as a function of p. When we do this, and use the same 
dimensionless factors as above, with the introduction of p = (ay/rnk) w for dimensionless 
momentum w, we have: 

dq{s) w{s) 

ds ~ ^l + 2aw{sY 
dw{s) ^ q{s) + q{sr) ^20) 
ds ^l + 2aw{sY 

V I — 

- (s - Sr) = V2cr\q{s) + q{sr)\. 

c 

The initial conditions are w{0) = 0, g(0) = 1, and we assume (for the purposes of computing 
Sr) that these values hold for all s < 0. It is clear that in the limit a — > 0, we recover 
the familiar instantaneous displacement, s = Sr, and the correct non-relativistic equations 
of motion. 

The dimensionless factor v/c gives us an additional knob - we can elect to set f <^ c 
while maintaining finite relativistic energies (the v/c — > limit functions, in the retardation 
condition, much like the a — > oo limit). The other limit, f = c is the one on which we 
will focus on the grounds that: 1. we have no physical model with which to fix v, and 2. 



using V < c allows multiple solutions to (16) (since the particles themselves can travel faster 
than V, allowing multiple points of causal contact). This could lead to interesting dynamics, 
however it is not clear how to weight the effects of the multiple solutions. Inspired again by 
electrodynamics, we will make the simplifying assumption that our field travels at the speed 
of light. 



We are ready to numerically solve the set (20 ) for arbitrary a. The method we use to solve 



the equations of motion in (20) is basically velocity Verlet |6], a finite difference approach 



where we use a fixed grid Sj = j As, and approximate the ODEs in (20 ) at these grid points 



For the retarded time condition in (20), we use a root-finding routine (bisection) adapted to 
the grid that will allow us to solve for Sr at Sj up to errors of size As/2. The scheme is easy 
to develop and we present it briefiy in the Appendix, so that the interested student could 
implement it. 



When the method (A2) is iterated, we generate approximations to the delay differential 



equation defined by (20), and the trajectories for a few values of a are shown in Figure 4 



What we see there is amplitude growth: as time goes on, it looks as if energy is being added 
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FIG. 4. The numerical trajectories approximating solutions to (20), shown for increasing value of 



a. In all cases, we use a step size, in (A2) of As = 0.01. 



to the system. Since we are using the relativistic second law, that additional energy makes 
the trajectory look more and more like a sawtooth. The result is quite striking, and given 
that it is a numerical solution, we would like to carry out some analytical tests to confirm 
the basic behavior. We will consider two checks, both appropriate for a small, where exact 
results are obtainable. Both of these apply to the relativistic well. 

The growth behavior of q{s) is evident even in the non-relativistic limit with a small, 
where its source is relatively easy to trace. The equation of motion becomes 

Cpq 



ds"" 



- [q{s) + q{Sr)] 
11 



(21) 



in our dimensionless variables - if a is very small, then s — ~ ^Jg is close to zero. We 
will approximate this small difference with a constant < e ^ 1 (a delay > 1 would occur 
only if we allowed the masses to move with speeds greater than c). We can then expand the 
equation of motion to get 

= - W) + q{s - e)] ^ -2g(s) + eg'(s), (22) 

and see that it has an "anti-damping" term eg'(s). The solution is 

(23) 



a 



a cos I - as \ — esm -as 



2 



a = Vo — e 



In this form it is clear that for a constant small delay the solution grows with time (at least 
for small s). 

As a further check, we can take a continuous trajectory, q{s) = cos(-\/2 s) (the solution to 
the non-relativistic problem with instantaneous forcing), and calculate both the equal-time 
spring force. Fa = — 2g(s), and the retarded time spring force, Fb = — [q{s) + q{sr)], by 
solving the retardation condition for this analytic function. The result is shown in Figure [5] 
where we use a = 0.1. What we see is that the retarded force is large (relative to Fa) as 
the particle nears equilibrium, so that the mass gets a significant kick near g = 0. Then the 
retarded force diminishes, and the particle continues with less deceleration than it would 
feel under the influence of Fa- The process repeats; each time the particle nears equilibrium 
Fb gets large compared to Fa and then eases off. The series of kicks, followed by decreased 
deceleration force, leads to an overall increase in amplitude. 

Physically, the growth must be caused by whatever is mediating the interaction - in E&M 
we have fields to transmit the effect of one particle on another, and those fields carry energy. 
We have no model for the corresponding mediating "field" in the retarded spring setting, so 
we have necessarily excluded its contribution to the total energy - we then have an "open" 
system, which doesn't conserve energy. 



VI. CONCLUSION 



The description of a relativistic oscillatory system is less universal than its non-relativistic 
counterpart. The immediate complication is the identification of ~kx as a component of 
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FIG. 5. The equal time spring force, Fa (in black), and the retarded time spring force, Fb (in 
gray) for a sinusoidal trajectory. 



the ordinary, three-force, or the Minkowski four-force that provides a starting point. If a 
fixed electrostatic source is generating a harmonic potential, then the force we use in the 
relativistic Newton's second law must be 

Ft = -kx (24) 



where A; is a constant (associated with the electrostatic sources) and x is the position of a 

dp 
d7 



particle in the lab. But we could also imagine —kx as appearing in ^ = —kx, leading to 



Fr = -kx^l-^. (25) 
We see that each is a reasonable candidate if we were to introduce a "Hooke's law force" in 



the relativistic form of Newton's second law (and there are many others defined by (10)). 
Both forces lead to the same limiting behavior in the low and high energy limits, while 
differing in the intermediate region. We can track these changes explicitly by examining the 
change in period as a function of a, the ratio of initial potential energy to rest energy. 

We introduced a time of flight delay in the separation x appearing in F^- in the sim- 
plest possible way. The resulting retarded force causes amplitude growth which seems non- 
physical. In a closed dynamical setting of particles generating and responding to fields, the 
force we have introduced would lead to energy loss by the field, but we have no way of 
performing the accounting that would balance the total energy of the system (lacking the 
field portion of the story). 
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Aside from developing some effective field with its own dynamics, it would be interesting 
to introduce damping in the retarded setting to try to tame the amplitude increase (that 
damping could come, for example, from radiation reaction if the particles are charged). 
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Appendix A: The Numerical Method 



Our goal is to discretize the dimensionless equations of motion (20), and develop a re 



cursion relation that allows us to move forward in time from the initial conditions. Let 
q{sj) = qj and w{sj) = Wj, then from Taylor expansion we know that 

dq{s) 



ds 



2 As ' 



(Al) 



and similarly for We proceed to define an update by replacing derivatives in (20) with 



finite approximations, obtaining 



Qj+i = (Ij-i + 



2 As Wi 



w 



[Sj - Sr) 



2 Ask [qj + q{sr)] 



(A2) 



l + 2aw] 



^\qj + q{sr 



Given starting values go and wq, we can use (A2) to propagate the solution forward in time. 



At each step, we must solve the retarded time condition by finding k such that 



[sj — k As) = \/2a \qj + q^l 



(A3) 



using the set of positions {qe}izl. In practice, we cannot achieve equality in (A3) using our 
fixed grid, so we take the k that comes closest to equality. That rounding introduces error 
As in the determination of the retarded time, and hence As in the evaluation of the force 
associated with the position at the retarded time. Then the update for Wj^i is accurate up 
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to errors that go like As^ per step (since the retarded evaluation q{sr) appears with a factor 
of As in front of it in (A2)). 
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